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STANDING ALFVEN WAVE CURRENT SYSTEM AT lo : ' 
VOYAGER 1 OBSERVATIONS - 


M. H. Acuna, F. M. Neubauer and N. F. Ness 


Abstract 


The enigmatic control- of the occurrence frequency of Jupiter’s 
decametric emissions by the satellite Io has been explained theoretically 
on the basis of its strong electrodynamic interaction with the co-rotating 
Jovian magnetosphere leading to field aligned currents connecting Io with 
the Jovian ionosphere. Direct measurements of the perturbation magnetic 
fields due to this current system were obtained by the GSFC magnetic field 
experiment on Voyager 1- on 5 March 1979' when it passed within 20,500 km 

south of Io. An interpretation in the framework of Alfven waves radiated 

t ' 6 

by Io leads to current estimates of 2.8 x 10 amps. A mass density of 7400 

3 

to 13600 proton mass units per cm is derived which compares very favorably 
yith independent observations of the torus composition characterized by 7-9 
proton mass units per electron for a local electron density of 1050 to 1500 

_*3 

cm . The power dissipated in the current system may be important for 
heating the Io heavy ion torus, inner magnetosphere, Jovian ionosphere and 
possibly the ionosphere or \ even the interior of Io. 


Introduction 


The Voyager 1 spacecraft encounter with the planet Jupiter occurred on 
5 March T979 with a closest approach distance .of 4.9 Rj .(1 flj = 71372 Km) 
at 1.2:04 UT. The encounter trajectory had been designed specifically to 
provide a passage through the magnetic flux tube which links the Galilean 
satellite Io with the Jovian atmosphere and ionosphere, at close range to 
the satellite. Distortions of the. flux tube associated with Io’s 
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interaction with the Jovian magnetosphere were not considered in the 

i 

trajectory optimization analysis due to the absence of definitive 
information regarding the magnitude of the expected perturbation. The 
selected trajectory was designed to provide a passage through Io's southern 
flux tube at a closest approach distance of 20500 km from the satellite. 

The twin Voyager magnetic field experiments (Behannon et al., 1977) 

consist of an array of four triaxial fluxgate magnetometers designed to 

measure magnetic fields over an extremely wide dynamic range, from the 

noise level of the low field ’magnetometers (LFM* s) , 0.006 nT RMS, to the 

6 

upper limit of- 2x10 nT per axis for" the high field magnetometers (HFM's). 
Two LFM’s are mounted on- a 13 meter boom and arranged in a dual 
magnetometer configuration which allows the continuous estimation and 
elimination of spacecraft, generated magnetic fields; the twc HFM's are 
-.mounted on .the boom cannister support structure and are designed to measure 
large planetary fields, up to 2x10^ nT per component. 

During the time interval associated with the anticipated Io flux tube 
passage, the GSFC magnetic field experiment detected significant, perturba- 
tions superimposed on the much larger background planetary field. We 
interpret these perturbations as associated with intense currents electro- 
dynamically induced in Io and/or its ionosphere as the Jovian magnetosphere 
sweeps past the satellite. 

. A preliminary report and analysis of quick look data obtained at the 
time of the encounter has been given by Ness et al. (1979). These pre- 
liminary analyses did not consider the ramifications of the observations 
reported by complementary investigations aboard the spacecraft. Of 
particular importance are the plasma inertial effects introduced by the 
presence of the Io plasma torus (Broad foot et al., 1979; Warwick et al., 
1979; Bridge et al . , 1979). It is the purpose of this paper to present a 
more detailed and expanded report of the magnetic field observations and 
the results of analyses carried out taking into account a more accurate 
physical description of the interaction between Io, its environment and the 
Jovian magnetosphere and using updated and much more accurate spacecraft 

trajectory and especially attitude information. 

$ 
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We show that the observations can be interpreted within the general 
framework of theoretical predictions based upon two oppositely directed 
currents in- each Jovian hemisphere flowing along the direction of Alfven 
wave group propagation, approximately parallel to Io's magnetic flux tube 
and generating a field perturbation whose general characteristics can be 
well represented by the field associated with a two dimensional dipole. 
This more accurate interpretation of the observations in terms of a 
non-linear, standing Alfven current system (Drell et al., 1965; Neubauer, 
1,979) yields estimates of the local ion mass density. These are in 
excellent agreement with the results derived from the observations by the 
ultraviolet experiment (Broadfoot et al., 1979), the planetary radio 
astronomy experiment (Warwick et al . , 1979) and the plasma science 
experiment (Bagenal et al., 1980; Belcher et al, 1980). In addition a 
number of interesting conclusions can be derived about the general na'ture 
of the electrodynamic interaction in the vicinity of Io. 

Earlier Studies 


The role of Io as a unipolar generator and the induction of large 
electrical currents flowing aligned to the magnetic field along the flux 
tube had been suggested on observational and theoretical grounds by 
Piddington (1967) and Piddington and Drake (1968). Further studies on the 
resulting current system, its possible closure geometry and on charged 
particles in the flux tube and their acceleration, have been presented by 
Goldreich and Lynden-Bell (1969), Gurnett (1972), Goertz (1973), Goertz and 
Deift (1973), Deift and Goertz (1973), Shawhan et al. (1975), Dessler and 
Hill (1975, 1979), Cloutier et v al. (1978), and Dessler and Chamberlain 
(1979). The strong theoretical interest in this unique phenomenon is 
directly associated with the observed correlation between Io’s. orbital 
phase, the Jovian Central Meridian Longitude (CML) and the probability of 
detection of Jovian radio emissions at decam etric wavelengths at the earth 
(e.g., see recent review by Carr and Desch, 1976 and references therein). 
These phenomena and the recent discoveries by Voyager 1 of the Io plasma 
torus (Broadfoot et al . , 1979) and of active volcanism on the surface of 
the satellite (Smith et al . , 1979) undoubtedly make Io one of the most 
interesting objects in the solar system. 
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A wid’el.y -accepted model-- for- the-- current system consists of -field- ■ 
.aligned currents emanating from- the -outer -face .of Io and flowing both 
southward and-, northward ' towards the Jovian ionosphere. The return- currents 
then> .flow- along the flux .tube on its .inner face. The closure of these 
currents at the extremes of the circuit requires the existence of, -regions 
of high conductivity transverse to the magnetic field. These were 
generally -associated, with the conductivity of the Jovian ionosphere, Io'.s 
ionosphere -and/or its* interior-. 

The description- of the interaction- in terms of a- non-linear standing 
Al-fven- current system- (Neubauer, 1 980). also, implies- the possibility that 
these waves- may- be . reflected at the Jovian- ionosphere or torus boundary. 
This would lead- to, a partial closure- of the- current -system in the near, 
vicinity of Io. In addition, in the- Alfven-wave system -description, the 
curttents are- not strictly field aligned and the separation between current- 
stwfeam lines and magnetic field lines can-become appreciable if followed 
a&i the'- way -down to the Jovian ionosphere (Neubauer, 1980). This -canrhav.e 
important implications for the interpretation of Jovian decametric radio 
emissions- as a function-»of CML and io-^phase since it further .enhances- the 
distortion of field lines away from -simple- meridional planes -associated 
Viph-the- non-dipolar- nature of the Jovian magnetic field (see., for example, 
Goldreich and Lynden-Bell, 1969). 

Observations and Analyses 

t 

, - / 

A projection- of the spacecraft encounter trajectory onto the Jovian 

magnetic equatorial plane in System- III (1965) coordinates is shown in' Fig. 
1. The position of Io is also indicated- in the figure for reference. The 
Io L-shell is crossed during the inbound leg at 0900 SCET while; the 
predicted Io flux tube traversal occurs on' the outbound leg at- 1502 SC.ET. 
One -interesting feature of this trajectory . is that the spacecraft remains 
relatively cl-ose to Io for 8 hrs,. thus potentially allowing- the detailed 
study of various phenomena (i.e., ion cyclotron waves) associated- with the 
Io torus. These studies will, be the subject of a future report and are not- 
included here. 



• The selected flux tube geometry relative to the spacecraft trajectory 
was calculated using the 0^ magnetic field model of Acuna and Ness (1976) 
and the model of Smith et al. (1975), neglecting any possible 
distortions of the flux tube due to the current system. The probability of 
intercepting the central part of the flux tube, as computed from the above 
models, was essentially identical for both, the only difference being the 
encounter time, with the 0^ model predicting a flux tube encounter ^ 1 
minute later than the model. The detailed trajectory near lo's closest 
approach is shown in Fig. 2 in a coordinate system centered at and moving 
With lo. In this analysis system, the Z-axis is parallel to Jupiter's 
rotation axis, the X-axis is oriented in the direction of corotational flow 
and the Y-axis points toward Jupiter. The closest approach point to the 
satellite occurs at 1513 SCET at a distance of 20500 km from lo's center 
and at a System 111(1965) longitude of 328°. 

To ensure the acquisition of maximum resolution data during the lo 
flux tube passage, the outboard LFM was commanded to operate in a fixed 
(non-automatic ranging) mode corresponding to a maximum measurement 
capability of 2400 nT per axis, close to the maximum anticipated for this 
time period from the analytical models of the field. This operational 
mode, in conjunction with the experiment 12-bit A to D converter and fast 
sampling .rate, yielded a data set with 0.5 nT and 60 msec resolution 
respectively. 

The magnetic field data obtained from 1424 SCET to 1550 SCET are shown 

in Figure 3, where 1.92 second averages of the three field components, B r , 

B' , B in a spherical coordinate system centered at Jupiter and aligned 
6 $ 

with the rotation axis are shown. It is important to note that the 
spacecraft performed a number of attitude maneuvers after 15:22'. SCET 
associated with optimum positioning of the scan platform for observing lo. 
The spacecraft attitude for some of these maneuvers has been reconstructed 
from the ‘magnetic field data itself due to .the lack of adequate high time 
resolution engineering data . 

During and adjacent to the time interval associated with the 
geometrical flux tube passage, strong, smooth perturbations of the 

9 
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■background, field, can be observed:. Although 1.92 -second averages are shown 
in the fig.ur.e., detailed analysis of the high time resolution data reveals 
the absence .of short period coherent perturbations in the measured field. 

A .power spectral analysis of thes.e data, up to' the instrument's Nyquist • 
frequency confirms the absence of significant maxima in the spectra, the 
power spectral density vs,, frequency .curves being characteristic of 
quantization noise associated with the ± 0.5 nT digitization uncertainty 
for frequencies greater than - 1 'Hertz. 

In order to study the characteristics of the perturbation, the. large 
background field must be subtracted from the measurements. Although 
analytical models of the .main magnetic field of Jupiter can -be utilized in 
principle to represent the field, this approach does not take into account 
the diamagnetic effects associated with Io's plasma torus and/or 
perturbation fields associated with the Jovian equatorial current sheet 
t which are significant - at this radial distance, considering the magnitude of 
lo's current system signature detected -by Vo.yager 1. 

Estimates of the components of the local background' magnetic field 

* / 

have been obtained by performing a quadratic least, squares regression 
analysis of the data excluding the interval where the principal Io induced 
perturbations are observed. This approach- yields a more accurate local 
■representation of the background field for the purpose of its subtraction 
fron the measurements. After this is carried out, the resulting perturba- 
tion field components AB R , AB g and clearly indicate that the observed 
perturbation is nearly transverse to the amb.ient magnetic field, as shown 
in Fig. 4. Since the spacecraft latitude is' close to the Jovian equatorial 
plane, the measured B Q component is approximately aligned with the ambient 
field direction. Note that the largest perturbations are observed in B R 

and B ppinting to the transverse nature of the perturbation field. This is 

<t> ' . , ' ' 

further illustrated in the top panel of Fig. ^ -where the angle bet.ween the 

background, and the perturbation field has been plotted showing, that it 

approaches 90° during this time. The computed perturbation vectors exhibit - 

a weak dependence upon the selected intervals for the calculation of the 

background field and this leads to the quoted- uncertainties in the 

estimation of parameters for the model described below. 
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At this point,' we can raise the question about the observed perturba- 
tion being due to an intrinsic Ionian magnetic field which is strong enough 
to be observed at the spacecraft location. Two simple dimensional 
arguments can be easily put forward to rule out this possibility: 1) the 
magnitude of the transverse perturbation vector shown in Fig. i| as a 
function of time varies with a characteristic scale length at its full- 
width- half-maximum points corresponding to a distance of 9000 km, as seen 
by an observer on Io, which is small compared with the closest approach 
distance of 20,500 km and 2) the perturbation is characterized by a change 
in ambient field direction rather than field magnitude. On the other hand, 
the possible existence of a small intrinsic Ionian field (Neubauer, 1978; 
Kivelson et al . , 1979) cannot be ruled out on the basis of these observa- 
tions. The existence of an intrinsic magnetic field of Io would physically 
affect the current system responsible for the observed perturbation in the 
immediate vicinity of the satellite (Neubauer, 1980) and would cause 
energetic particles to be selectively deflected as a function of their 
•energy. This latter effect has been reported by Krimigis et al . (1979) 
although no definitive statements have been made regarding the possible 
existence of an Ionian magnetic field. 

Qualitatively, the extremely smooth variation of the pertubation field 
together with its simple appearance suggest that Voyager 1 did not pass 
through the main current carrying region but flew near it. As reported by 
Ness et al. (1979), the observations can be represented to a good 
approximation by the magnetic field generated by two oppositely directed 
line currents flowing approximately parallel to the average magnetic field 
direction. 

Because of the infinite nunber of possibilities available to explain 
the given magnetic field observations along the spacecraft trajectory in 
terms of a spatial distribution of currents, detailed interpretations of 
the observational results have to be guided by appropriate physical ideas 
about the nature of the plasma currents. Further insight into the geometry 
of the current carrying region can be gained by considering the inertia of 
the ambient plasma (Io torus) which will cause the currents to deviate from 
flowing strictly parallel to the magnetic field. This problem has been 


7 



addressed recently by Neubauer- CI 980) where a full non-linear solution' is 
given for the ‘system of standing Alfven waves generated by. the interaction 
of Id with- the co-rotating, magnetosphere. The current tube in this case 
forms an angle 6^ with the ambient magnetic field and is determined by the 
local Alfven characteristic direction, as illustrated in Figure 5. Ibis 
description of the interaction- suggests that the observations should be 
analyzed in a coordinate system aligned with the direction of Alfven 
characteristics C* (see Fig-. 5)‘. Since the magnitude of depends on the 
local Alfven Mach number and this is not known _a priori , it must be 
assumed’. Thus, the orientation of the analysis coordinate system with' 
respect to the local magnetic field i's determined by' the assumed A-lfven 
Mach number, which yields a local characteristic direction 1 by 


1 + = 


% 
T&i 


+ m a * 


•d) 


where % is the 1 ' average background field vector, ie the local Alfven tech 
number = V ^/V^) and $ the unit vector in the $ direction of the 
spherical coordinate system centered at Jupiter. V pe ^ and are the speed 
of Io relative to the m'agneto spheric plasma and the Alfven speed' 
respectively. Equation (1) reflects the fact that the $ direction 
represents the direction of corotational flow. 


We choose to orient the 2-axis of our Io-centered analysis coordinate 
system antiparallel to the direction of Alfven characteristics' and 

define. the X-Z plane to contain the direction of corotational flow. The 
Y-ax-is then completes the right-handed system and is roughly directed 
radially towards. Jupiter . The magnetic field perturbation data and 
spacecraft positional information are then rotated into the analysis 
coordinate system. Note that because equation (1) the unperturbed 
background magnetic field is not aligned with one of the coordinate 
system axe’s. Therefore the perturbation magnetic field has three non-zero 
components in this representation. 

The coordinate system chosen i-s most suitable for model fitting. 
Except fo.r the immediate vicinity o^f the - source, Io, the physical 
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quantities aB, b$, etc. do not depend on the 2-coordinate and the magnetic 

field components aB , aB can be represented by a scalar potential y(x,y) 

x y 

if the current component along the Z-direction, j_ = 0. These currents j 
(x,y) feed the currents in Io r s inmediate vicinity or are fed by them. AB 
is then given by the requirement of constant magnnetic field magnitude. The 
component AB z (x,y) is due to the relatively small currents forming closed 
loops in the X,Y-plane. Strictly speaking, the general current streamlines 
are therefore helices with small pitch angles. 

The lowest order terms in the expansion for the potential ¥ of the 
field generated by a system of currents j (with net current = 0) flowing 
parallel and anti-parallel to the 2-axis in our coordinate system is a 
2-dimensional dipole whose potential, 'K is given by 

♦d) . ° ‘ d - S) to 

2n C [X| - |D| 2 ) 

and + 

AB = - y 7i|i (2) 


where corresponds to the perturbation field, X is the spacecraft 
position vector X(x,y)., m is the dipole moment and D the dipole location in 
the X-Y plane. 

The observations can then be fitted to the model in the analysis 
coordinate frame in a least squares sense. Given the Alfven Mach number 
which defines the orientation of the coordinate system, we determine the 
parameters m and B by means of =an iterative "walking dipole" algorithm 
which minimizes the root-mean square sum of the residuals. Since the 
coordinate system is centered at Io and the current must flow approximately 
parallel to the Z-axis to connect with the satellite, the dipole location 
obtained must be as close to the origin as possible (note however that the 
curvature ‘of the field lines is ignored although the fly-by distance of 
20,500 km is not insignificant). If a dipole location reasonably close to 
the origin is not obtained, the assumed Alfven Mach number is incremented 
and the fitting process repeated for each Alfven Mach number chosen. 
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Ibis- process was' carried oat for.Alfven Mach numbers in the range of 
.01 < *25. • The best fit to the observations in terms o.f the 

connection -of the current system" to Jo is obtained- for an Alfven- Mach 
number" of 0.15' ± .0-1. The- perturbation field vectors along- the spacecraft, 
trajectory in the Alfvenic coordinate system are shown in Fig.. 6 for = 
0.15. Also, shown are the orientation and location of' the 2-d'imensional- 
ddp.ole: jnonent m obtained from the fit" to the observations.- The orientation 
of the d-ipole, approximately parallel to the direction of corot'atibnal 
flow, is in. excellent agreement with- theoretical expectations, that- is, a 
downward current flowing along the outer face of the current tube while the- 
upwards return current flows- along the- inner face. The’ Y-coordinste- of the 
dipole location is not exactly’ zero but exhibits a small negative value 
which agrees -reasonably well in magnitude and direction with the- fact that 
the small spatial variations of the global background magnetic field of 
Jupiter, i.e. its curvature, have^been : ignored in this model. Other 
mechanisms which could account for small displacements of the center of' the - 
current system include the effects of finite frequency corrections" to the 
MHD- equations , an asymmetry of the interaction region near lo, etc. The- 
circle; in the figure represents the projection of Io on the X-Y plane and 
the dipole moment m obtained is 8.1 ± .2 x 10^ amperes km. THi-s is 
consistent with the earlier estimates (Ness et al . , 1 97'9 ) but' at their 
lower limits. If it is assumed 'that the currents flow on the periphery of 
a cylinder- of Io's diameter (3680 km) according to a sin- 4> distribution, as ' 
indicated in the figure, the derived" current flowing along- the Alfvenic 
tube is 2.8 ±'.1 x 10^ amperes. 

The Voyager 1 magnetic field observations constitute the first in situ 
experimental verification of the existence of a-current system induced by 
the interaction of Io with the magnetosphere. However, since the 
observations are restricted to the the spacecraft trajectory and its single 
traversal of the fl-ux tube region, no- unique conclusions can b'e drawn from 
these observations regarding the location of the current closure regions, 
at Io or in the Jovian ionosphere, their geometry, or any possible 
dependence of the intensity of the currents upon the relative location of 
the flux tube with respect to Jovian longitude. 
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It is important to note that the quantity derived from the 

•observations is the- dipole moment m and not the current. A few remarks on 

the uniqueness problem are in order at this point. Although the above 

choice of a current distribution is a reasonable one for physical models 

where the current closes in the interior of Io or its adjacent ionosphere, 

there is an infinite number of current distributions which produce the same 

field as an idealized two-dimensional dipole along the fly-by trajectory in 

the X,Y-plane. This is apart from the ambiguity due to observational 

uncertainties and the fast decrease of higher multipoles as a function of 

distance from the origin. In addition, it cannot be concluded that the 

current carrying region is remote from the spacecraft trajectory from the 

fact that 3 ' two-dimensional dipole is the lowest order term, i.e. the far 

field term-, in the expansion of an arbitrary distribution of current j z 

with net current zero. For example, a cylinder of diameter 8700 km which 

just touches the Voyager trajectory and carries a sinusoidal current 

distribution' with total current = 1.2 x 10 amperes, yields exactly the 

same' perturbation field as the current distribution above. The smallest 

total current consistent with the observations is obtained by a sheet 

current located just inside the fly-by trajectory as seen from Io. Its 

value is given by I . = 6 x 10 amperes. Both examples illustrate the 

min. 

non-uniqueness problem.. The final justification for a current distribution 
must ccrae from auxiliary observations and physical modeling. 

Figure 7 illustrates the excellent agreement between the observations 
(■48 second averages) and the results derived frcm the model. The small 
discrepancy which is evident after 1520 SCET is not completely understood 
at the present time. The inclusion of terms of higher order than the 
dipole in the expansion for the scalar potential il> does not appear to be 
justified to explain this discrepancy since the model fits the observations 
extremely well elsewhere prior to 1520 SCET. The possibility that these 
observations reflect seme temporal effects cannot be ruled out at the 
present time and further analysis,, taking into account correlative 
observations by other experiments on, the spacecraft, is required. There 
may also remain some uncorrected spacecraft attitude errors. 

• Discussion and Summary 
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The results presented above, based upon more accurate trajectory and 
attitude dpta, support the interpretation of the Voyager 1 magnetic field 
observations neat Io (Ness et al., 1979). They show that the spacecraft 
did not traverse the current carrying region itself but passed approxi- 
mately 4350 km from its center. The resuits given here take into account 
the high plasma densities implied by the presence of the' Io' torus. Strong 
support, for the interpretation of the observations in terms of a standing 
Alfven current system is given by the fact that the Alfven Mach number 
required to connect the currents to Io utilizing a 2-dimensional dipole 
model is 0.15 ± .0.1,.= This value is in excellent agreement with that 
derived from the results reported by the planetary radio astronomy and 
ultraviolet experiments. Warwick et al.0979) and Birmingham et al. (1980, 
this issue) have estimated the electron density in the torus at this time 
as 1050-1500 el/ cm 3 . Taking into account the Alfven Mach number derived 
from the magnetic field observations, = .15, we can estimate the ion 
mass composition in the torus as 7-9 proton masses per electron, in 
excellent agreement with the results reported by Broadfoot et al. (1979) 
for the plasma composition S^, O^^) . Fran the value of the Alfven 

Mach number we can also derive the external Alfvenic conductance Z A = 

1/(u V.) (Neubauer, 1980). In addition, the Alfven wave theory predicts a 
simple relation between the magnetic field and velocity perturbations: AV 
- Al/(v 0 p) 1/2 . For example, the maximum magnetic field perturbation of 70 
nT yields |A$| = km/ sec. It will be interesting to test the 

predictions of the theory with the plasma ion observations during the Io 
encounter. The predicted velocity perturbation vectors' have been indicated 
in Figures 6 and 7. Preliminary comparisons with plasma observations ( 
Belcher et al . , 1980) indicate .good agreement with this picture. 


Neubauer (1980) has derived a maximum magnetic moment which can be 
generated by the current system as 


m 

max 


2 * R Io E o M A 


„ 2 . 1/2 

*o (1 + M A } 


(3) 


12 



where we have assumed 8=0, corresponding to transverse flow (see 

Neubauer, 1980 for derivation), B q is the ambient magnetic field and 2 R^ o 

is the current system diameter. For B 0 = 1900 nT, = 0.15, Rj q = "WO km 

we obtain m max = 4.77 x 10^ amperes km, which is smaller than the value 

derived from the observations. As pointed out by Neubauer (1980), this 

would imply that the current carrying region may have larger dimensions 

than the assumed lo diameter and that the ionosphere of lo may play a more 

significant role than its interior in providing a closure path for the 

9 

currents. To obtain the derived moment of 8.1 x 10 amperes-km, the 
'effective radius' of lo in equation (3) has to be increased to 
approximately 2400 km. This would imply that the current carrying region 
is located 560 km above lo' s surface, this value should be compared with 
the results obtained by Iiliore et al., 1974 for the extent (750 km) and 
scale height (200 km) of the daytime Ionian ionosphere. It is interesting 
to note that volcanic ejecta can reach altitudes of 270-300 km during an 
eruption (Smith et al., 1979) and thus could contribute a significant 
amount of ionization at these heights through collisional processes with 
the charged particle population in the radiation belts or the torus plasma 
An alternate explanation for the large magnetic moment may be the 
non-negligible effect of reflected Alfven waves in the case of very 
efficient short-circuiting of the electric field by lo (or non-stationary 
processes). For this case to occur the short-circuited electric field 
. inside the current tube must be below an upper limit given by the 
condition that the round trip travel time of the Alfven waves to the 
reflecting boundary (torus or Jovian ionosphere) equals the time it takes 
for the plasma -to traverse an lo diameter. In a torus of diameter 2R . we 
obtain 


‘W M A * 


W 


which represents a generalization of equation (29) of Neubauer (1980) and 
is also valid for Alfven waves in the region of strong short-circuiting. 
Here E o is the corotational electric field (see also Goertz, 1980). Using 
s 0.15 we obtain E^/E^ = 0.17. Should the reflected waves be 

insignificant we must have E i /E o > 0.17. Since the average value of the 
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•density and therefore in the torus is probably somewhat larger .than 
sbet.ween -Io and Voyager, the transition electric field .may be even smaller 
■’than this value.;. 

The implications of the derived current intensity of 2/8 x 10^ amperes 

for the electrodynamic interaction of -Io with its environment are 'extremely 

interesting, particularly the .power ’dissipated in the northern and southern 

portions of the current system. In deriving this quantity we have to take 

into account that the various closed current filaments are subjected to 

differing voltage drops. Therefore an average voltage drop 'has to be 

multiplied by the total current. We 'Obtain P = 2jmx"6 I. -Using "£ = 

^rel x ^o’‘^ an< ^ B o = ^900 nT a val’ue of P = 1.8, x 10 watts is obtained 

which is independent of. the diameter of the -current carrying cylinder and 

therefore less model dependent than the total current. This value of P, as 

pointed by Ness et al . , 1979', is comparable to the energy dissipated "in 

lo’s interior by tidal forces (Peale and Cassen, 1979) and the energy 

radiated by the plasma torus at ultraviolet wavelengths (Bro’adfoot et al . , 

1979). For the total power to be available for heating Io or its atmosphere 

or ionosphere it. would be necessary that Io and its environment provide the 

largest resistance in the total current circuit. This would be in 

contradiction to the observation that the magnetic, moment m is close to its 

maximum value or even exceeds it. Essentially total short-circuiting of the 

electric field in Io’s vicinity is required’ for m (observed) - in to. be 

max 

attained and therefore negligible resistance in Io and its neighborhood. We 

therefore conclude. that most of the power is dissipated outside Io and its ' 

immediate environment , i.e. the torus, the' inner magnetosphere or the 

Jovian ionosphere. We estimate the maximum power dissipation in id’s 

atmosphere or possibly its interior due to Joule heating to be-P/10 = 2 x 
11 

10 watts. 

The currents flowing through the Io system will also give rise to an 

accelerating torque due to the rotation of the Jovian magnetic field, as 

* 

first discussed by Goldreich and Lynden-Bell, (1969). The value of this 

torque, calculated using (T = 2 r .m.B Y with r T = the distance of Io 

T6 io o io 

from Jupiter, is 1.3 x 10 Newton-meters which is insignificant compared 

35 p * 

to Io 1 s orbital angular manentum of* 6.5 x 10 kg.m /sec. Thus, 



conclusions reached by Goldreich and Lynden-Bell (1969) regarding the 
characteristic doubling time for Io's angular momentum of 10 years stand 
unchanged. (Note that these authors used values of 3500 nT and 1.1 x 10 5 
amperes for the magnetic field and current intensity, respectively), lie 
also mention that only a fraction of the torque is generally applied to the 
solid body of Io.' For example the portion of the current near Io due to the 
pick-up of newly generated ions by the magnetic field (Goertz, 1980) will 
not contribute to a torque on the solid body if these ions experience no 
further collisions with Io’ s surface or atmospheric neutrals after 
generation. Currents in a collison dominated region of the neutral 
atmosphere will contribute to the torque via viscous forces. 

Prior to the discovery of the Io plasma- torus, it was generally 

accepted that the magnitude of the currents flowing in the flux tube was 

limited primarily by the Jov.ian ionospheric conductivity (Cloutier et al., 

1.978; Dessler and' Hill, 1979). The magnetic field models developed after 

the Pioneer 10 and 11 encounters with Jupiter (Smith et al . , 1975; Acuna 

and Ness, 1976) indicated that the field was asymmetrical and that the 

surface field intensity exhibited broad minima at certain longitudes where 

enhanced particle precipitation would be expected (Roederer et el., 1977; 

Dessler and Hill, 1979). This led Dessler and Hill (1979) to postulate 

that the Jovian ionospheric conductivity at the foot of the Io flux tube 

would be greatly increased at these longitudes due- to the enhanced particle- 

participation, and this, in turn, would lead to an enhancement of the Io 

flux "tube currents. This ’’active hemisphere” or "magnetic anomaly" model 

predicts that the Io currents should increase to over 10 5 amperes in the 

range of System III (1965) longitudes of 200-260°, compared to a value of'' 

./'10 5 amperes over all other longitudes (Cloutier, 1978; Dessler and Hill, 

1979). The results deduced here frctn the Voyager 1 observations at a 

System III longitude of 310° < X < 335°. corresponding to the Io flux tube 

crossing interval, show that this current is 2.8 x 10 5 amps rather than the 
5 

expected -e 10 amps. Since the magnetic anomaly model does hot take into 
account the existence of the plasma torus, the role played by it in 
providing local closure to the currents must be fully evaluated. The 
observations reported here suggest that Io acts as a generator of strong 
Alfven waves with its associated currents, magnetic field and velocity 
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perturbations. In. contrast to the classical models by Piddington U96^),- 
Goldreich and Lynd_en-Bel’l (1969) etc. the disturbances issuing from Io 
during, the. Voyager-1 encounter wer'e most likely determined by the local 
torus properties i .while the interaction with the torus boundary or Jovian 
ionosphere was small' (a more detailed discussion of this point is contained 
in Neubauer , 19,80). ,We also conclude that most of the power dissipation in 
this current systfem occurred in the Jovian ionosphere, inner, .magnetosphere 
or torus but not inside Io or its itrmediate vicinity. 

From the picture presented above it is clear that, the coupling between 
Io and the Jovian ionosphere is not. as strong nor as conceptually simple as 
envisioned prior to the Voyager 1 observations and plasma tor.us discovery. 
This raises new and interesting questions regarding the mechanism of Io 1 s 
control of deeametric radio emissions. Although it could be concluded that 
-the lead angie provided by the Alfven wings (8.5°) - is a reasonable 
explanation for the observed Io-phase asymmetry in the probability of 
detecting deeametric emissions, little is known about the geometry of the 
postulated reflection of the Alfven wave at the Jovian ionosphere. 

The large increase in Alfven velocity outside of the plasma torus 
significantly reduces the magnitude of- this angle leading to a possible 
overlap of the downgoing and reflected waves. This effect would further 
complicate the reflection geometry at the foot of the flux tube where 
presumably the emissions originate. 
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Figure . Captions - 


Figure.'. 1., Voyager .1 trajectory and position of Io -in -System III longitude 
Cl 965) and radial, distance (in _ units of Jov.ian ,Radi-i)... on 5 ;Marchrl979.- 

Figune 2. Relative position of Voyager 1 and Io during, Io current system 
studies . 

Figure 3; GSFC magnetic . field observations in. spherical .coordinates. The - 
three orthogonal components' of the field are- .shown as is the 
magnitude. Data' points are 1. 92-, second averages' of 60 msec.-, vector 
samples . 

Figure Perturbation magnetic field associated, with the Io- current 

system, as detected by -Voyager 1. In addition to the three. orthogonal" 
components in spherical coordinates - the- angle between the background 
field and perturbation field has been computed in the upper- panel. 

Figure 5. Schematic diagram illustrating the, electrodynamic interaction of" 
Io generating Alfven waves in the ambient plasma torus.- The specific 
Mac'h numbers chosen are for illustration purposes- only. and do not' refer 
to the results obtained by Voyager 1 except that the flow is 
sub-Alfvenic and probably, not supersonic.. 

Figure 6. Perturbation field due to Io current system as observed by 
Voyager 1 along its trajectory and the position of the equivalent 
dipole moment , m, of' the current system. Also shown is the°current 
density vector along a circle with Io’ s radius. The- predicted scale for 
velocity perturbations is also- indicated. See text. 

Figure 7a. Comparison of magnitudes of magnetic field perturbation due to 
Io current system as observed' (+) and predicted (continuous curve) from 
dipole line currents model. 

Figure 7b. Same as Figure 7a. but a different fitting interval for the 
background magnetic field has-been used. In this case the value of 

' t r ./i 

obtained is 0.16. 
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